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TELECOMUNICAÇÕES 


O planeamento do complexo sistema 
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cutado pelo Departamento Técnico 
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laboração com a Comissão de Elec- 
trificação da C. P. 
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Estes rolos são apenas uma pequena fracção da nossa gama de fitas 
adesivas tesa. 

Temos muitas qualidades diferentes: incolor e a cores, em celofane 
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aplicação. 
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ROCK MECHANICS IN PORTUGAL “* 


RESUMO 


E apresentada uma aintese du actividule que tem sido 
desenvolvida em Portugal no dominio da Mecânica das Ro- 
chas, Ieferem-se quer vs estudos em laboratório quer «in 
situ», saltentando-se as contribuições para o progresso dus 


conhecimentos e das técnicas e vs resultados mais relevantes 


by MANUEL ROCHA 


Director, Laboratório Nacional 
de Engenharia Civil, Lisboa 


SPNOPSIS 


A synthesis of the activities in the field of Roch Mecha- 
nica in Portugal is presented, The article describes labora- 
tory and «in situ» studies, specially those contributing to 
the progress of knowledge and the improvement of techni- 


ques, anel the major resuls 80 far obtained, 


que têm sido obtidos, 


1. INTRODUCTION 


In Portugal, activities in the field of Rock Mechanics have been almost exclusively concen- 
trated at the Laboratório Nacional de Engenharia Civil (LNEC). These activities began about twenty 
years ago, in connexion with the foundations of Castelo do Bode dam, and have been progressively 
developed since then. 

Studies in the field of Rock Mechanics have been stimulated above all by the problems posed 
by the foundations of concrete dams |!'. Nevertheless opportunities have also arisen for studying 
problems concerned with underground power stations, tunnels and slopes. 

At present the staff active in Rock Mechanics studies at the Laboratório comprises 20 full- 
-time employees, 10 of which are university graduates. Additionally, other sections of the Labora- 
tório such as those dealing with instrument design and construction, observation of the behaviour 
of structures, and documentation, help Rock Mechanics studies. 

The policy followed in Rock Mechanics research is the same adopted by the LNEC as a 
technological research institution. Thus: i) a constant attention has been paid to the identification 
of problems with a real practical interest, which led us to give a special emphasis tho the in situ 
study of rock masses and to the observation of the behaviour of structures ; ii) it is attempted to 
obtain a balanced knowledge, i.e. the available means are not concentrated on particular aspects of 
a problem at the cost of neglecting others, perhaps as relevant ; iii) the improvement of knowledge 
and the accuracy of predictions are adjusted to the real needs of the problem to be solved. 


(*) Apresentação da actividade desenvolvida em Portugal no domínio da Mecânica das Rochas, feita no 
I Congresso da Sociedade Internacional de Mecânica das Rochas, que teve lugar em Lisboa, em Setembro de 1966. 
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I am going to describe the activities so far developed and those in progress. Due to the short 
time available, no more than a general picture of the aspects deemed more important can be presented. 


2. QUALITY INDICES 


Given the marked heterogeneity usual in rock masses, in particular near the surface due to 
weathering, we have been seeking for some time to establish quality indices which, by means of a 
quick and economic test, can characterize the properties of a rock mass, e.g. its deformability or its 
shear strength. This is possible only if, in each case, a correlation, to be obtained by means of in 
situ tests, exists between the index and the property to be assessed. 

In the case of granular rock masses little influenced by jointing, as frequently occur in grani- 
tes when weathering is much advanced, it has been possible to obtain correlations — sufficiently 
tight to be of practical interest — between the mechanical properties of the rock mass and the water 
absorption by specimens of the rock :?]. The correlation between the modulus of deformability 
determined in situ and the water absorption, considered as a quality index, is presented in fig. 1 
for the granite rock of the foundation of Alto Rabagão dam, a structure about 90 m high recently 
completed in Northern Portugal [!]. The quality index was determined at many points at different 
depths as the excavation proceeded, until values corresponding to the deformability modulus assu- 
med in the design were obtained at the different foundation zones. Thus by means of a quality 
index it was possible to decide more objectively than usually, the required foundation depth. 
Notice that on the basis of model tests an extremely low deformability modulus of about 10 000 kg/cm* 
(140 000 psi) was accepted in the design for the decomposed granite forming the upper portion 
of the left bank. This value was confirmed by observing the displacements of the dam during the 
first filling of the reservoir. 

As it was to be expected, we have been unable to find significant correlations between the 
quality indices and the mechanical properties of rock masses when the behaviour of these is appre- 
ciably influenced by jointing and other macrofractures. In the fig. 2 the abcissas are the values of the 
quality index and the ordinates the values of the modulus of deformability determined in situ, 
during the exploration of the foundation rock mass of the future Vilarinho das Furnas dam. As can 
be seen there is no correlation, 
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Another quality index that we recently began to use, which is dealt with in a paper presented 
at this Congress [*], is based on the micropetrographic analysis of thin sections of the rock. The 
micropetrographic index is defined as the ratio, for a given volume of the rock, of the volume of 
unaltered minerals to the volume of altered minerals, fissures and voids. A correlation obtained 
between the moduli of elasticity of core samples and the micropetrographic indices is presented in 
fig. 2. As seen, this correlation is comparatively satisfactory. As this index takes into account not 
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only the voids, as does the preceding index, but also the altered minerals, it is to be expected that 
it yields tighter correlations with the mechanical properties of the rock. Anyway the interpretation 
of the mechanical behaviour of rocks in their mineralogic composition and of the alteration of the 
minerals seems very promising. 

Quality indices are especially helpful in the case of rocks from which specimens for mechanical 
tests cannot be obtained, due to their alteration or their structure, v.g. schistosity. This is parti- 
cularly important in civil engineering undertakings, in which it is usually attempted for economic 
reasons to reduce excavations to a minimum, so that comparatively superficial formations have to 
be considered in areas where, due to weathering, undisturbed samples are difficult to obtain. 

When specimens for mechanical tests in laboratory can be obtained, it is better to select for 
quality index a mechanical property, such as the modulus of elasticity, the compressive strength or 
the shearing strength, since the correlation is tighter between the values of this property and the 
mechanical properties of the rock mass. 

* The problems of alteration and alterability of rocks recently began to be studied and we 
intend to investigate in particular the influence of the state of stress. The problem of petrification 
of soils, particularly of lateritic soils, has been the object of some studies [*]. 


3. DEFORMABILITY OF ROCK MASSES 


I will now describe the activities in connexion with the three basic problems posed by the 
characterization of a rock mass, that is the investigation of its deformability, its strength and its 
state Of stress. In accordance with the general policy mentioned, in which we try to investigate as 
fully as possible the real behaviour of rock masses, a special attention has been paid to the study 
of these three problems by means of in situ tests [1] [*]. Both our and others” experience shows 
that the deformability and the strength of rock masses cannot be satisfactorily determined, as a 
rule, by means of laboratory tests. 

We usually determine the deformability of rock masses by means of jackloading tests in adits, 
in which an area of one square meter is loaded as shown in fig. 3. So far, more than 500 tests 
were carried out with this setup on formations of very different natures, so that a comparatively 
extensive knowledge on the deformability of rock masses was obtained. In special cases of rock 
masses in which heterogeneity and cracking are not very marked, jackloading tests are sometimes 
carried out on smaller areas, for instance circles 30 cm in diameter, with a view to preliminary 
estimates of the deformability. 

Creep is systematically investigated with the same setup by means of tests with a suitable 
duration, sometimes of more than a week. The results of in situ creep tests under stresses of about 
50 kg/cm? on areas of 1 m* [!] are presented in fig. 4. In abscissas are plotted the moduli of defor- 
mability corresponding to a load applied during 30 minutes, designated by Eso, and in ordinates the 


values of the relative decrease of the modulus of deformability, that is EM — Bom in which Es is 


Eso 

the modulus of deformability for an infinite time, obtained from the tests by assuming a certain 
creep law. The tests performed so far on rock masses of very different natures, grouted or even in 
the natural condition, show that the relative decrease of the modulus of deformability is less than 
15 º/o in the majority of cases and only in one case did it reach 25"/o. That is, creep is not rele- 
vant under stresses of the order of magnitude of those applied by the structures on their founda- 
tions. Notice that this has been observed even in very altered rock masses with moduli of elasticity 
of less than 5000 kg/cm”, which is rather unexpected. 

The tests for determination of the deformability are very often carried out on the rock mass 
in its natural condition and repeated on the same surfaces, after the rock mass has been grouted [!]. 
In fig. 5, the values of the modulus of deformability of the rock mass before grouting, En, are 
plotted in ordinates and the values of the same modulus at the same points after grouting, Es, are 
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E 
plotted in abscissas. Many of the points lie near the straight-line *- = 1, which means that defor- 


n 
mability remained practically unchanged with grouting. In some cases, however, the modulus of 


E 
deformability increased three times, as some points lie near the straight-line —* = 3, As is well- 


n 
-known, the effect of grouting depends on the nature of the rock masses, and very frequently their 
deformability is not improved by grouting. 
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The procedure used for determining the deformability of a rock mass is open to two criticisms: 
i) the tested zone of the rock mass is disturbed by the excavation; ii) the volume involved is not 
yet representative in many cases, although the loaded area amounts to 1 mí. 

In order to overcome these difficulties a new method, presented in fig. 6, is being developed 


for determining the deformability of rock masses. A slot or a series of contiguous slots are cut in 
the rock by means of a diamond disk with a diameter d, equal to for instance 1 m, and a width 


A flat jack is installed in each slot and the same pressure p is simultaneously applied in all the 
jacks. The increase in the distance between the faces of the slots due to the applied pressure p is 
measured by means of instruments inserted in the jacks. Fig. 7 shows the slot-cutting machine fitted 
with a disk 1 m in diameter. Although not yet applied in the field, this method seems an important 
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improvement as, on one hand, it practically eliminates disturbances in the rock mass and, on the 
other, it enables large volumes of a rock mass to be quickly and economically tested. 

We believe nevertheless that the jackloading tests will remain in use for the cases in which 
it is just the deformability of the rock mass, taking into account the disturbance due to excavation, 
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that is to be determined. Even so, we intend no. m a 

to improve it by measuring the displace- Hitieás O NTE DISK DRIVE 
ments of the rock mass at different depths, Ni fai O INLET ES | Pina RANA 
instead of only at the loaded surface, so as dt 
to obtain the deformability in the successive 

layers, where the disturbance decreases in HAIN 


function of the distance to the surface. It DRNE. | 
| 


- Ol HYDRAULI HOLDING TUBE 
« MOTOR | 1] 


DIAMOND DISK 


will thus be possible to predict, for instance, 

Jhe deformation of a foundation by integra- 

ting the deformations of different layers. E 
Another method in use for determi- + a ma de 

ning the deformability of rock masses, in 

special for designing tunnel linings, is the Fig. 7 

well-known method consisting in determining 

the deformation of a section of an adit under 

water pressure [*]. The setup used is presented in fig. 8. The deformations of the adit are measured 

along different directions by means of vibrating wire strain meters. The change in length of 3 dia- 

meters in function of the applied pressure is presented in fig. 9. In spite of a previous grouting, 

the rock mass presents a marked anisotropy with a larger deformability along the vertical, as is 

usual in altered rock masses. 


ROCK 


The methods so far mentioned for 
determining the deformability require 
that the point of the rock mass to be 
tested be made accessible by means of 
an adit which is a severe limitation, in 
special for points at great depth and rock 
masses with an intense water percolation, 
particularly at points under river beds. 


é Ri % CRESSURE CAGE That is why a dilatometer was recently 
2 WATER INLET 8 VIBRATING WIRE METER developed for determining the deforma- 
: e uno e rolar E bility along boreholes. As this dilatomne- 
5 WATER METER O INVAR RODS ter is described in a paper presented at 
6 WATER OUTLET 1º ADE this Congress [6], we will merely state 

Ná that it was developed for use in boreholes 


with the usual diameter in site investiga- 
tions. Thus, as shown in Fig. 10, the 


instrument has a diameter of 74mm RADIAL DEFORMATIONS 


with a view to enabling its use in e. OF THE GALLERY 
boreholes 76 mm in diameter. The n E 

deformations due to the pressure o 

applied on the wall of the borehole x e 

are determined very accurately by e 

means of differential transformers uu 

along four diameters at 45, so as to S 


investigate the anisotropy of the rock 
mass. The influence of cracking of the DIAMETER VARIATION (mm) O 008 mm 
rock mass on the deformation of the 
diameters was studied theoretically. 
Curves of thepressure versus deformations of two diameters, obtained in one test in which succes- 
sive loading-unloading cycles were applied, are presented in Fig. 11. As can be seen, the curves 
are entirely analogous to those usually obtained in jackloading tests. 


Fig 9 
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The dilatometer has the drawback that 
only a very small volume of the rock mass 
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TRANSFORMER 
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is tested. Nevertheless, given the possibility siuo o riNG SPRING LIQUID UNDER RUBBER 
of making many tests, almost amounting to |. ó di PE ti 
a continuous determination of deformability = H WO = = E era) ba 
along the boreholes, it should be possible, “ (dg Ade A E 
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obtained to a statistical processing, by means io CUM pIPE FOR THE LIQUID | 
of which the average deformability of the -755 sida o ii 
different zones of the rock mass could be de. femndass Fig 10 
termined with the accuracy currently required. 
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systematically carried out for stu- 
dying the deformability of rock 
masses. In fact, besides contribu- 
ting to a general characterization 
of the component rocks, these tests 
often enable a satisfactory study 
of the influence of certain parame- 
ters on the deformability of the 
rock mass proper. This is possible 
in the case of parameters with a 
slight influence on the contribu- 
tion of joints and other fractu- 
res to the deformation of the rock 
mass. Such can be the case, for 
instance, with the creep and the 
influence of the humidity on de- 
formability. 


Different laboratory studies have been carried out on the anisotropy of rocks. One of these, 
presented at the Congress |7] concerns the law of variation of the modulus of elasticity and of 
the ultimate compressive strength in function of the direction of the applied stress. The tests on 
granites so far carried out show that the distribution of these magnitudes is ellipsoidal. In fig. 12 are 
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Fig. 13 


shown the points representing the moduli of elasticity and the ultimate strengths determined in 9 
prisms taken from a single block along the directions indicated in the figure. As shown, the expe- 
rimental points differ but little from the solid line curves lying on the ellipsoids. According to the 
tests, the axes of the modulus of elasticity and of the ultimate strength ellipsoids do not coincide 
at each point, and the anisotropy of the modulus of elasticity has been systematically more marked 
than that of the ultimate strength. 

As for anisotropic rocks, I will mention further the study of the rheologic behaviour ofs chists 
and of the laws governing the change of deformation with the direction in an orthotropic body 
made up of alternate thin layers of two elastic materials. Both studies were presented as papers 
at this Congress |*| [']. 

We have also studied creep phenomena in rocks by means of laboratory, mainly triaxial, 
tests. A system of 9 triaxial cells for quick tests amd creep tests of rocks under confining 
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pressures up to 1000 kg/cm” is presented in Fig. 13. Specimens with diameters up to 7.5 cm can 
be tested, and longitudinal and transverse strains are measured by instruments applied on the spe- 
cimens inside the cells. These triaxial cells were developed for studying creep phenomena in rocks 


under the high stresses occuring in the neighbourhood of deep tunnels, 
4. MODELS AND FIELD OBSERVATION 


In order to take into account the influence of the deformability of the foundations on the 
behaviour of structures, particularly om their state of stress, model tests have very often to be 
used. Such is the case with concrete dams |!| ||. Fig. 14 shows the model of an arch dam 
whose foundation presents a system of holes whose diameters and distance were varied from 
zone to zone in order to reproduce the deformabilities in the rock mass as determined by in situ 


rock tests. In fig. 15 is shown a model of Tang-E Soleyman dam in which the sedimentary layers 
forming the foundation were reproduced |!). 


Fig. 15 


The study of more than 200 models of concrete dams of different types, in the majority 
of which the influence of the deformability of the foundations was investigated,enabled the conclu- 


sions presented in fig. 16 to be drawn. The figure expresses, in qualitative terms, the influence of a 


LITTLE LITTLE VERY 
IMPORTANT IMPORTANT IMPORTANT IMPORTANT 


change in the deformability of the fcundation on the state of stress in the dam, in respect to the 


case of a foundation with the same deformability as the concrete of the dam. The ratio —- of the 


Ec 
modulus of elasticity of the rock mass to that of the concrete can vary, even in extensive zones of 


the foundation, from high values down to 1/4 without an important influence on the state of stress; 
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contrary to some beliefs, it is only for values of É <— that the influence of the deformability 


Ec 4 
Er 1 E; 1 
becomes of importance, and considerably more so for — << — ; for values of — < — the 
c 8 e 


conception of the dam and the state of stress in it are entirely governed by the deformability of 
the foundation. These conclusions are of course of a general order and purely qualitative, as the 
influence of the deformability of foundations depends on numerous parameters, variable from case 
to case. 


Fig. 17 


The conclusions presented show that, save for foundations with very low moduli of elasticity, 
it is sufficient to determine an order of magnitude of the deformability of foundation rock masses. 
This means that accurate methods are not required; what matters is to obtain significant values, 

It is wo th noting that faults, re- 
levant joints and other geologic features 
are usually reproduced in our models. A 
model of an arch dam in which several 
faults were reproduced is shown in fig. 17. 
Notice that, provided a failure by sliding 
along a fault does not occur, the state of 
stress in a concrete dam of any type is, 
as a rule, only very slightly influenced 
by the presence of faults even if their 
width is considerable. 

Still regarding the influence of dis- 
placements in the foundation, a technique 
was developed at the LNEC for applying 
in models the action of earthquakes, 
which are repoduced as random vibra 
tions, as they actually are [**]. The setup 
used is presented in fig. 18 which show 
s the electromagnetic generator of random vibrations, with the required distribution of spectral den- 
sity of acceleration, coupled to the model of a building supported on 4 columns. The instruments 


Fig. 18 
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for measuring stresses, displacements and other magnitudes at points of the model can also be seen. 
As said we pay great attention to the observation of the behaviour of trusctures ["]. 
In fig. 19 is shown a sketch of a setup for measuring deformations along adits in rock mas- 
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Fig. 21 


ses, on basis of several meters. The variation of the distance between two points is determined by 
readings of reference marks on an invar-tape by means of microscopes connected to the points. It is 
thus possible to obtain an accuracy of 5>< 10º in the determination of strains. Two types of strain 
meters with long bases, whicg are applied inside boreholes in rock masses, are shown in fig. 20; 
these strain meters have been used in observations inside rock foundations of dams. 


Fig. 22 
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For measuring the vertical displacements at different depths we use the setup presented at 
the left part of fig. 21; at the right is shown an optical plumb that we intend to use for measuring 
the two components of the horizontal displacement of points along a drill hole. 


5. SHEAR STRENGTH 


1 will now describe the activities connected with the second basic problem of the characteri- 
zation of rock masses: how to determine their strength. As is known, it is the shear strength as a 
rule that is of more interest in the equilibrium of rock masses. 

The characterization of shear strength has been based mainly on in situ tests [1] [14]. Fig. 22 
presents the setup we use, in which blocks with a 70 > 70 cm cross section are tested. In order to 
increase as much as possible the influence of the heterogeneity and of the cracking of the rock 
mass on its shear strength, the distance d is given the highest value for which tensile bending stres- 
ses do not occur. In fig. 23 is showa a machine recently developed for preparing the lateral faces 
of test blocks by means of a diamond disk. It is thus possible to reduce the disturbance in the rock 
mass, which with the usual techniques of preparation of the blocks can entirely conceal the results 
of the shear tests. 

More than 250 in situ tests so far carried out with the described setup, not only on rock 
masses but also on joints and faults, form an invaluable body of experience. 

Fig. 24 shows the values of the coeffi- 

is cient of friction and of the cohesion in an al- 

tered granitic rock mass in function of the 
water absorption of the rock, taken as qua- 
lity index. The rock mass in question presents 
highly variable degrees of alteration, as the 
water absorption ranges from 3 to 30º/, 
from zone to zone. As can be seen, whereas 
the coefficient of friction decreases from 1.9 
to 0.72, i.e. a reduction of about 2.5 times, 
when the quality index increases from 3 to 
30, cohesion changes from 13 to 0,3 kg cm, 
i.e. is reduced as much as 40 times. J men- 
tion this fact because it is very general, i.e. 
the variation of the coefficient of friction 
along a rock mass is usually much less mar- 
ked than the variation of the cohesion. This 
is very important, because it shows why, in 
the analysis of the stability of rock masses, 


Fig. 23 


the coefficient of friction can have a much lower 
factor of Safety than the cohesion. 

As for the factors of safety to be adopted 
in each case, we had occasion to study the problem 
having concluded that as a rule a concept of sa- 
fety based on a reduction of the values of the me- 
chanical properties of the rock mass is preferable 
to the usual concept, according to which the 
applied loads are assumed to increase until failure 
occurs |). 

From the hundreds of in situ shear tests car- 
ried out we conclude that as a rule the angle of 
friction, even in very altered rocks, exceeds 50º and 
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even 55', i.e. is fairly, high whereas cohesion, save in sound rocks, is very low [!]. This tact, toge- 
ther with the high factor of safety required for cohesion, leads us usually to neglect cohesion in the 
rock masses and, on the other hand, to use sometimes angles of friction appreciably higher than 
those commonly adopted. 

As for the joints, in situ and laboratory tests have yielded rather low angles of friction even in 
the case of joints without clay infillings. In 16 in situ tests of joints occurring in schists and gra- 
nites, 11 values lower than 35º were obtained. The coefficient of friction correspording to 35" is 0.7, 
that is the value usually adopted in analyses of stability. As rock masses are usually divided by 
numerous joints, the results just mentioned clearly show the major role played by the interlocking 


Fig. 23 


of the blocks of which rock masses are made up, as otherwise the majority of foundations would 
not be stable. This indicates how important would be a theory enabling this interlocking effect 
to be calculated. 

As said, shear tests are also carried out in laboratory. The machine used, in which specimens 
with cross sections up to 20><20 cm can be tested, is shown in fig. 25. We have recently started 
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Fig. 27 


systematic comparisons of the values of the angle of friction in the joints as determined in situ and 
in laboratory tests. So far this comparison has shown a fairly satisfactory agreement between both 
values. The determination through laboratory tests of the angle of friction of the joints in rock 
masses themselves would be of great practical interest. 


Fig. 28 
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6. DETERMINATION OF THE STATE 
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4 o X - Ao Something will be said now about 


= ç ns 7 flag 
xk,| 5412 au x h ZÁ q] =—s the activities in connexion with the third 


ts E - fm! problem posed by the characterization of 
24 = rock mass: namely, the determination 
“CMATIERTEr a 1 
' K Xo of its state of stress. 
“1 In tunnels and underground power 


stations the state of stress has been de- 
termined in the walls by a stress relief 
method consisting in taking cores out by 
BOREHOLE ' > means of drilling or manual tools and in 
Bed em, 4 ABADE measuring the deformations with resis- 
tance strain gauges [!º]. The results obtai- 
ned have indicated the presence of tec- 

tonic stresses. 
STRESS-RELIEF Due to the difficulties in the appli- 
Y-BOREHOLE cation of this method and to its limita- 
+ tions, a new technique was recently de- 
veloped for applying the well-known 
Fig. 30 stress restoring method, which is descri- 
bed in a paper presented at this Con- 


CYLINDER WITH 
STRAIN GAUGES 


gress ["']. Its new feature is the fact that the slot is cut with a dia- 
mond disk, the segment obtained being shown in fig. 26. The slot 
cutting machine developed is presented in fig. 27; it is fitted with 
a disk 60 cm in diameter operated by an electric motor. The seg” 
ment-shaped flat jacks to be installed in the slot for restoring the 
initial state of stress are presented in fig. 28. The jacks fit exactly 
in the slot without need of mortar or other material. After the test 
the flat jack is removed without difficulty. A series of flat jacks are 
used, by means of which stresses can be determined as the depth 
of the slot increases. 

The method developed has the following advantages: i) it 
reduces disturbances in the rock mass to a minimum ; ii) it is very 
rapid and inexpensive; iii) it makes it possible to detect any ano- 
malous deformations while the slot is being cut; and iv) it enables 
the simultaneous determination of the modulus of deformability 
of the rock mass. 

Determination of stresses in the walls of tunnels has a 
great interest as regards the state of stress near the tunnels 
themselves. But when the state of stress is to be obtained at a 
point inside a rock mass, such determination has two serious 
well-known disadvantages, namely, the need to consider stress- 
“concentration factors around the tunnel, and above all the distur 
bances which the excavation of the tunnel creates in the rock 
mass. That is why a new method is being developed enabling the 
complete determination of the state of stress along a borehole. 
Considering a borehole in a rock mass along a direction OZ 
(Fig. 29), the method enables the three normal stresses 7x, 7y, Tz, 
and the three shear stresses, Tyz, 77x, Txy, to be determined. The Fig. 31 
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Fig. 32 


method essentially consists in (Fig. 30): à) drilling a hole with a small diameter, for instance 
3.7 cm; ii) cementing on the wall of the hole a plastic cylinder in which are embedded 9 electric 
resistance strain gauges; iii) relieving the stresses in the rock mass around the cylinder by over- 
coring with for instance a 7.5 cm diameter. The cylinder is glued to the walls of the hole yb means 
of araldite which, previously placed at the bottom of the hole, is extruded/when the cylinder contai- 
ning the strain gauges is fitted into the hole. Fig. 31 shows at the left the strain gauges mounted in 
a transparent pipe inside which the plastic cylinder is moulded ; at the right is seen one cylinder 
with the gauges embedded. A plastic cylinder inside a rock core is shown in fig. 32. In order to 
assess the accuracy of the method it was applied to prisms of rock maintained under a constant 
compression in a testing machine. The results obtained with a limestone prism, under a stress of 
100 kg/cm?, are presented in fig. 33. It can be seen that the values of the strains measured along 
the 9 directions indicated in the figure do agree with the expected strains, and that the components 
of the state of stress, in respect to the system of axis Oxyz, calculated from the measured strains 
are in a still better agreement with the imposed stresses. 


O = 100 kg/cm? 


STRAINS (107º) STRESSES (kg/cm?) 
Expected Measured Calculated | Imposed 
sn oe At a si ; 
Ej == Es — 195 — 183 Cs 50 50 
E = €9 26 3 330 C, | 031 O 
X €3 63 71 C; 53 | 50 
6 =6 | — 102 — 105 | Tx | 0043) O 
Es 491 48 8 To | 0021 | 0 


ER 200 som | Tre 424 | 50 


Fig. 33 


7. FUTURE ACTIVITIES 


Finally I would like to present our plans for the future. We intend to go on developing our 
activities, particularly in connexion with specific problems of tunnels and underground power sta- 
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tions. Furthermore we believe it is advisable to increase the proportion of our means invested in 
basic research, as the existing experience and knowledge ensure the efficiency ot such an activity. 

Given the present status of Rock Mechanics, we will go on paying a special attention to what I 
shall call the inquiry on the real behaviour of rock masses, i.e. their in situ characterization, and 
the observation of the behaviour of the structures built on or inside them. In parallel with this 
activity, the studies in laboratory of the properties of the rocks will go on, having always in mind 
their relationship with the properties of the rock masses. As another item of our future activities I 
would mention a more intensive use of models not only for predicting the behaviour of rock masses 
in definite cases, ie. in connexion with given structures, but also for investigating phenomena 
in general. 

The theorization of the knowledge obtained in both field and laboratory investigations will 
be carried forward as far as possible. 
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O MÍNIMO DE PRODUÇÃO DE ENTROPIA NO ESCOAMENTO 
UNIFORME DUM FLUIDO NEWTONIANO 


RESUMO 


Demonstra-se que no escoamento laminar e uniforme 
dum fluído Newtoniano a condição de mínimo de produção 
de entropia é equivalente às equações de Navier-Stokes, 

Usando a formulação variacionul das equações do mo- 
vimento que decorre da demonstração anterior mostra-se 
como podem obter-se soluções analiticas exactas ou aproxi- 
madas para o escoamento laminar, uniforme, em condutas 
de secção constante. Como cxemplo trata-se o caso duma 
secção rectangular para o qual se deduzem as soluções 
exactas e se comparam as soluções aproximadas, 

Também para a secção rectangular se representa gráfi- 
camente o factor de resistência em Junção da relação entre 
os lados e se mostra como a utilização do factor de resis- 
tência deduzido para tubos pode conduzir a erros muito 
upreciáveis, 


1 — INTRODUÇÃO 


por  J. DELGADO DOMINGOS 


Prof. 1. S. T. 
Do NEEM do L. A.C. 


SYNOPSIS 


The author demonstrates that for a Newtonian fluid in 
laminar, uniform flow, the condition of minimum entropy 
production is equivalent to the Navier-Stokes equations, 

Useng the varcational formulatron, resulting from the 
former demonstration, it is shown how exact or aproximate 
analytical solutions for laminar, unitorm flow in conduits 
with constant section, can be obtained, As an example, the 
exact solutions for a revtangular section are calculated 
and compared with the aproximates ones. 

Also for a rectangular section, the friction factor is 
plotted as a function of the sides ratio, and vt 18 shown 

how the use of the friction factor deduced for pipes can 
introduce important errors in the results. 


Foi demonstrado por Prigogine [1] e posteriormente generalizado por De Groot |2] na termo- 
dinâmica irreversível linear que para sistemas contínuos em equilíbrio mecânico (!) e possuindo coe- 
ficientes fenomenológicos constantes a produção de entropia é mínima em regime estacionário. 

A nosso conhecimento, porém, não se tratou ainda, pelo menos explicitamente, o caso de sis- 


temas abertos com v fO (?) e tendo em conta a produção da entropia devida ao aparecimento de 
tensões tangenciais de viscosidade (*). Este aspecto tem para as aplicações na engenharia um parti- 
cular interesse, não apenas pelo seu conteúdo físico intrínseco mas também como meio natural de 
introdução das técnicas do cálculo variacional na solução exacta ou aproximada de problemas de 
grande complexidade analítica. O presente trabalho insere-se nesta via, na medida em que pretende 
mostrar como algumas técnicas simples do cálculo das variações permitem a solução analítica aproxi- 
mada de certo tipo particular de problemas em mecânica dos fluidos e o sentido físico intrínseco 
que possui a sua formulação. 


(!) No sentido restrito em que a velocidade dum elemento do contínuo não varia no tempo, acompanhando-o 


no seu movimento fm -— o) o 
Dt 


(?) No caso de sistemas com mais de um componente onde se verifica difusão o mesmo pode permutar massa 
sendo v=-0 [2]. 


(*) Note-se que pode ser simultâneamente rot v=0 e div v£0 o que dará origem à produção de entropia 
se não é nulo o segundo coeficiente de viscosidade. Todavia neste caso não existem tensões tangenciais. 
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2 — MODELO FÍSICO 


2.1 — Na sequência da exposição admitiremos que: 


— O sistema é constituído por um fluido Newtoniano incompressível, de viscosidade 
constante. 

— O sistema evolui a temperatura constante e em regime laminar. 

— A fluição verifica-se num canal de eixo rectilineo e secção recta (A) constante. 

— O eixo oz coincide com o eixo do canal e designando por u a única componente 


- . 4 “e d u 
não nula da velocidade, só nos ocupamos da região em que o = 0. 
z 
— As únicas forças aplicadas por unidade de massa são as que derivam do potencial 
gravítico. 


2.2 — Nas condições do modelo, a produção de entropia por unidade de volume e de tempo 
será [3]: 


ds du N? du NV? (') 
no [GAS] 
dt | o dy 
sendo “ — viscosidade dinâmica 
s — a entropia produzida por unidade de volume devido ao trabalho irreversível das 


forças tangenciais 
t — tempo 


e o aumento de entropia por unidade de tempo no elemento de volume Adz será 


= [ q dAdz 
dt A 
donde 


ds . du ) du | 
2 D = dA = = — — dA 
ai dt dz alii ds 


2.3 — Da conservação da massa resulta também que, sendo 2 a massa específica 


3) G = | ou dA 
“ À 
e Ge D são constantes para todas as secções rectas do escoamento. 


3 — MÍNIMO DE PRODUÇÃO DE ENTROPIA 


3.1 — Se considerarmos G como um constrangimento imposto, da família de superfícies que 
verificam 3) haverá uma que torna 2) extremo se existir uma funcional u (x, y) tal que 


( 2 / 2” 
Jal odox dy 


| 2 | 2 
() O termo [e + 7a — q designa-se por vezes como função de dissipação. 
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O sã ã .ÃÕÃÃcÃá rã 


sujeito a 


5) f em dA 
A 


e obedecendo no contorno de A às condições impostas para u que caracterizam o problema físico 


6) f(u) = 0 


3.2 — Demonstra-se em cálçulo das variações [4], que 4) e 5) são equivalentes à determinação 
do extremo de 


ROC ETR 


sendo À uma constante a determinar posteriormente. 
As condições necessárias para a estacionaridade de 7) obtém-se das equações de Euler 4] e 
reduzem-se neste caso a 
ou ou 


8 + ee 24— — À p=o0 
) - Sa aê g 


Da família de superfícies de parâmetro À que satisfazem 8) e 6) é solução do problema a que veri- 
fica 5) a qual permite deste modo a determinação da constante arbitrária introduzida. 


3.3— O movimento do fluido rege-se pelas equações de Navier-Stokes que nas condições do 
modelo se reduzem a 


du |, du dp 
9 —— — "| TES qu ———s e 
pias CR) add 


sendo p a pressão hidrostática e g. a componente segundo oz da força por unidade de massa. 

As equações de Navier-Stokes, juntamente com as condições aos limites f(u) =0 determinam 
u univocamente (*). Desta unicidade resulta que a condição 5) é equivalente a impor um certo 
valor para o segundo membro de 9) do que resulta que ), parâmetro da família de superfícies que 
satisfazem 8) e 6) seja determinado por 5) como 


10) À == 2 (= se + 5: ) 
z 


3.4 — De 8), 9) e 10) conclui-se que a existência dum extremo da produção da entropia sujeito 
à condição de em cada secção se escoar o mesmo caudal mássico é equivalente à exigência de se 
verificarem as equações de Navier Stokes sujeitas sobre o contorno às mesmas condições. 

A condição necessária de extremo é também suficiente e poder-se-ia demonstrar, considerando 
a segunda variação de D, que esse extremo é um mínimo. 

Nas condições do modelo pode pois considerar-se como um teorema que: 


— O perfil das velocidades em escoamento uniforme é o que torna mínima a produção de entropia para 
um caudal imposto, 


(') Se as condições sobre o contorno envolvessem apenas as suas derivadas a solução seria indeterminada a 
menos numa constante aditiva. Todavia, mesmo neste caso se manteriam as conclusões que se seguem, as quais se 
preferiu obter do modo que se apresenta, embora sacrificando o rigor formal, porque mais simples e fisicamente mais 
intuitivo. 
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Dada a natureza do problema variacional, é imediata [4] a demonstração de que a formulação 
anterior é equivalente (º) a 


— O perfil das velocidades em escoamento uniforme é o que torna máximo o caudal para um valor 
imposto da produção da entropia. 


3.5 — O regime estacionário corresponde pois ao de menor dissipação de energia mecânica 
compatível com os constrangimentos impostos ao sistema. Atendendo a que 


o= fole (eo) oe 


é intrinsecamente positivo ou nulo, a ausência de dissipação corresponderá a =0 ou a u==cons- 
tante =c. 

No primeiro caso trata-se dum fluido perfeito. 

No segundo caso seria necessário que também no contorno u=c e o fluido comportava-se 
como um sólido indeformável (7). 


4 — APLICAÇÕES 
Tornar mínimo 
p Ju? 2 4 
mo ca=[ (e) +05) ]taso 
A dx dy a 
Sendo 
“e 


12) B=— — 
Es P Bz 


é, como se viu, equivalente à solução das equações de Navier Stokes sujeito à mesma condição sobre 
o contorno 


13) f(u) =0 


A vantagem da formulação apresentada não consiste porém, à parte o seu sentido físico, 
na redução do problema à solução das equações de Navier Stokes mas sim no tratamento directo 
de 11) por alguns dos métodos conhecidos do cálculo das variações sobretudo quando se pretendem 
soluções analíticas aproximadas (º). 

Um desses métodos, o método de Ritz [4] consiste em propor como solução uma função 


n 
14) F (x, Y, Cl ... Cn) sm Ci Vi (x; y) 


em que as 7; (x,y) são linearmente independentes e verificam 13) e determinar os c; de tal modo 
que I, seja mínimo. 

As funções qy à parte verificaram (13) e serem linearmente independentes são em princípio 
arbitrárias. 

Substituindo (14) em (11) resultará 


15) je É S a) S - 
fue (Bait + pin de Ro 


(*) O problema clássico, que deu o nome a este capítulo do cálculo variacional — problemas isoperimétricos — 
foi o da determinação da figura plana que possui a maior área para o menor perímetro, ou a que para um dado perí- 
metro possui a maior área; em ambos os casos trata-se do círculo. 

(') Se considerarmos este como um fluido de » = co, de ser D -£ O resultaria ser também u=c. 

(*) A obtenção destas poder-se-ia fazer sem a prévia formulação variacional, como seria por exempo o método 
de Galerkin. Todavia, conhecida a expressão variacional correspondente à equação dada não apenas se obtém um novo 
sentido físico para o fenómeno como em certos casos a solução é mais simples de obter. 
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ou 


16) [= eye + 4 fi ci=f (c) 
Í 
Lj=1..n 
onde 
dv o doi od 
A dx dx dy dy 
17) 


fi=2B f v;dA 
A 


como 1 é apenas função dos c; a condição de estacionaridade exprime-se por 


18) O em id Dad 
d ci j 


Lj=1..n 


Sistema de n equações lineares a n incógnitas que determinam os ci. 
Das relações anteriores pode verificar-se ainda que 


1 i 
Fmédio = BA Bic; 


Pode demonstrar-se que se os 7 constituirem um sistema completo [4] F=Zc;y con- 
verge para a solução exacta quando n aumenta (eventualmente quando n— c9). 

Do ponto de vista das aplicações, a ilação é que da forma escolhida para q; dependerá a 
necessidade dum maior ou menor número de termos para uma dada aproximação. Na escolha de ? 
pode ser-se guiado pelo sentido físico do problema, por analogia com situações análogas, por 
intuição ... 


4.1 — Como primeiro exemplo de aplicação, considere-se o escoamento entre dois planos à 

distância 2a segundo x, e infinitas segundo yz, cuja solução exacta é sobejamente conhecida. 
. —- . ER 

O problema trata-se a uma dimensão, e usando o método proposto consiste em tornar minimo 


19) =[" | ral + 2 Bu | dx 


Como condição ao limite consideramos a ausência de escorregamento sobre a superfície 


20) u=o0 para x=+a 


O anulamento da função proposta para x ==“ a obtém-se simplesmente considerando como 
factor comum a 9; o termo (a? — x?), 
Deste modo 


21) u=(a2-x) [utcx +... +cn x"] 


satisfaz as condições exigidas de independência linear e veriflcação de 20). 
Consideremos apenas a primeira aproximação «q 0 c=c=0. 
Obtém-se, substituindo em 19) e integrando 


| = E a* [u ci? + Boi | 
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e de 


I 
— =0=22u+B 
des 
B 
resulta = — — 
2 
donde finalmente, atendendo à definição de B 
1 dP 
22) = (— a tre) (2-0) 
24 dz 
Se tivessemos considerado mais termos na expressão 21) além do primeiro, obteríamos 
cg=...=€=0, o que resulta do facto de 22) ser uma solução exacta, como pode verificar-se. 
É também 
1 B 
Um = — — a? nine 
3 p 


4.2 — Como segundo exemplo considere-se um canal de secção rectangular e lados 2a e 2b, 
respectivamente. 


Como condição sobre o contorno será, como habitualmente : 


23) u=0 para = “e à 
y=+b 


Estas condições poderão ser satisfeitas tomando como factor comum das funções 7;, 
(xº— a?) (y?— b?). 

A simetria de u em relação aos eixos sugere que as funções 7; em produto com (y? — b?) 
sejam formas quadráticas. 

Proporemos 


u=(-a)(y'-b) [utoxºtoy +... +enx"y"] 


4.2.1 Considerando apenas a primeira aproximação (º), obteremos sucessivamente 


au = Eco a* b' [a? + b?] 
45 


Oy == — ab B 
9 
oq spp pbjup É PbB-o 
dci 45 + 
5 B 1 
A O A RE 
8 u (a? + b?) 
donde 
5 B 
24) “= — — E = = 
8 REESA didi 


(') Como imediatamente se conclui das expressões gerais 16) e 17) a parte mais trabalhosa do método é a 
determinação dos ai;j e de 3. Deste modo, se após a aproximação com n termos, desejarmos considerá-la com n -;-1, 
todos os coeficientes já calculadas se mantêm. 
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2 ki 
25) sucos, 
18 uu at+hb! 


4.2.3 — Em segunda aproximação consideremos apenas o caso em que a=b. 
Dada a simetria relativamente às diagonais será c) == c3 e virá 


= (a? — x?) (a? — y?) [eu +ca (x? + y?)] 


Obtendo-se (1º) 


c4 == — 0,29219 Es 
ua? 
c2 == — 0,059228 EA 
H aí 
donde : 
B 
26) u=— a (a? —- x2) (aº— y?) [0,29219 + 0,059228(x? + y?)] 
B 
27) Um == — 0,1404 a? — 
H 


4 — MÉTODO DE KANTOROVICH 


4,3 — Em vez do método de Rietz pode usar-se, em muitos casos com vantagem, o método de 
Kantorovich ou das equações diferenciais ordinárias, 

Neste método propõe-se a solução combinando funções conhecidas de uma das variáveis com 
funções incógnitas da outra. Reduzindo o integral duplo a um integral simples por integração na 
variável em que as funções são conhecidas e estabelecendo as condições de estacionaridade do integral 
simples obtém-se um sistema de equações diferenciais ordinárias que determinarão as funções 
incógnitas. 

Exemplificando o método para o caso em estudo e usando apenas a primeira aproximação 
proporemos 


28) u= (by?) £ (x) 

substituindo em 11), virá 
a b 

2) pe f de | elbyye pay e)+2B (by? E) dy 
— A “ -— 

e efectuando a integração em y 


a 16 8 8 
30) = [ [e ee [+ Bb] a 
dh - 15 3 3 


(1) Um problema inteiramente análogo no domínio da Elasticidade é apresentado em [6], pág. 283, relativa- 
mente à torção plana. Se for a função das tensões, e M o momento de torsão, demonstram-se fácilmente as seguintes 
relações 
B B M 


21G6 ? u 4AGO 


que são sempre válidas se P e u obedecem sobre o contorno a P=u=o0. 
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A condição de estacionaridade de 30) exprime-se pela Equação de Euler correspondente (!!) : 


a [4 € |-nEs! E- E b'B=o 


dx | 15 3 
ou 
dt 5, 5 Bs 
dx?  2b? 4 pb? 
Equação diferencial ordinária de coeficientes constantes cujo integral geral é 
Vi x —/5 * ' 
b 2 b 1 
f=a e. +c e —— — 
2 
As constantes de integração determinam-se de f == 0 para x == + a obtendo-se finalmente 
cs lo chVS E) 
=D 0 lg MEO 
E cos h + E) 
a 2 D/. 
donde 
T e (4/55) | 
31) u=—— — (by?) 1— ——— > A 
* E cos h (y/ Es 
EL 2 


doe hd io À tenda A 


4.4 — A solução exacta pode obter-se tanto pelo método de Rietz como pelo de Kantorovich, 
desde que se proponha para u uma expressão sob a forma duma série infinita de funções que for- 
mam um sistema completo no domínio. 

Usando o método de Rietz e propondo 


OO a] 


u= Cmn Cos 
m=1,3,5...n =1,4,5 2 a 2b 


X 
(4) No caso geral de da F [x, & (0), fe (x) ..., En o), Fi (o)... En (x)] dx 
o 
as condições de estacionaridade, Equações de Euler, são [4] : 
A Eis E ) CAS 
dx id Fi / d É, E 


o 


E (+) E 


dx dn dê 


que determinam as n funções incógnitas. As constantes de integração obtêm-se das condições aos limites. 
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que ôbviamente satisfazem o u==0 para x= ta y=+b obtém-se [5]: 


m-+-n sd 
So co a . 
— 1) * = ) 
33) u= > > cas Mcash, rea cos ——— cos + 
n=1,3,5... m=1,3,5... * gu mn (m+E n?) 2a 2b 
onde E om E 
b 
OQ OO 
256 a B 
34) e Um =—— mo o 
= = = mn (m'+Fn') u 
Usando o método de Kantovorich ter-se-ia partido de 
2] n* 
35) ue o cos * f(x) 
n==1,3,5 2b 
e obter-se-ia (12) o 
n-—! n=x 
' oo Z cos h “ab 
36) u=— Ana > SE) q ag Em cos 2") 
Ro UV n=185.. nº | cai RUE 2b 
E 2b 
e atendendo a que 
16 bº Ss (-1) dE n 
s nf 
- a * cos y 
nº p= 1,35. nº 2b 


E. . é ed he . - 
é o desenvolvimento em série de Fourior de — (b*— y”) no intervalo (—b, b), a expressão 36) assume 
2 


a forma mais simples: 


37) u=— — EE st) = 
u 2 


Por outro lado 


38) Um = — — 


e — 


4.5 — Comparação de soluções. 
4 5.1 — Para apreciação do valor relativo das várias soluções aproximadas faz-se seguidamente 


a comparação numérica do valor médio da velocidade e do seu máximo, usando as várias expressões 
apresentadas. 


(7) As expressões que se seguem correspondem às apresentadas em [6] p. 276 muito embora obtidas partindo 
da condição de que 33) seja solução da equação às derivadas parciais que exprime a estacionaridade de 11). 
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Valores de Kk= ——— 


es b? 
HH 
F = a/b 1 2 3 oo 
- 
Solução exacta 0,1402 | 0,2286 0,26533 0,3333 
Método de Kantorovich 
ê 0,1396 0,2283 0,2630 0,3333 
à áp. 
Método de Rietz 
á 0,1388 | 0,2222 0,2500 0,2777 
1.º ap. | 
Método de Rietz | 
. 0,1404 
2. ap. | 
(*) Supõe-se b fixo 
(0,0 
Valores de — 
— Bb*/u 
E = a/b 1 | 2 e 
Solução exacta 0,2938 0,1138 0,5000 
Método de Kantorovich 
0,3026 | 0,1144 0,5000 
1.º ap. | 
Método de Rietz | 
0,3125 0,5000 0,625 


1.º ap. | 


Método de Rietz | 
2.* ap. | | 


Como se verifica, usando o método de Kantorovich o caso mais desfavorável corresponde 
a F=1.0O erro em um é da ordem de 0,5'/, mas em umax da ordem de + 3 “A. ParaF >2 0 erro 
cometido é despresável, pois a aproximação converge rápidamente para a solução exacta. 

Para a primeira aproximação do método de Rietz o caso mais favorável é para F=1. O erro 
em Um é da ordem de 1 º/,) e em umax da ordem de 6 "/o. Para F>1 o erro cresce rápidamente até 
um máximo de 16,7 º/. É todavia apreciável a melhoria que se obtém com a segunda aproximação 
para F=1. O zero em um é da ordem de 0,15 "/o e em umax de 0,3 *f. 
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Como se verifica também, a 1.º aproximação de ambos os métodos comete um erro por 
defeito em um e por excesso em umax. Na segunda aproximação do método de Rietz verifica-se 
o inverso. 


4.5.2 — Dos processos apresentados infere-se que, na primeira aproximação, o método de Kan- 
torovich é preferível. Para secções com F = 1 a sua vantagem é incontroversa. Deve acentuar-se 
porém que uma conclusão de carácter geral não é possível pois o erro cometido depende essencial- 
mente da forma das funções escolhidas para solução. Todavia, na medida em que no método 
de Kantorovich a função incógnita corrige em certa medida o erro cometido com a função proposta, 
infere-se que seja o método recomendado sempre que não se possua à priori uma noção correcta 
da forma que a solução final deve assumir, se a determinação de função incógnita não se revelar 
de dificuldade excessiva. 

Como é óbvio, o valor da solução aproximada depende da finalidade que se pretende. Se esta 
for o conhecimento do valor médio de u é em geral da maior utilidade mesmo a primeira aproxi- 
mação. É o caso por exemplo de apenas de desejar o cálculo dos factores da resistência. 

Se é o conhecimento de u em todo o domínio a finalidade da solução aproximada, os desvios 
podem ser importantes nalgumas zonas. No caso presente devem esperar-se os maiores desvios junto 
aos cantos. 

O conhecimento de u e não apenas o seu valor médio ou máximo tem particular importância 
em transmissão do calor porque o seu conhecimento prévio é necessário para a determinação 
do campo das temperaturas. Todavia, mesmo neste caso, a utilização de soluções analíticas aproxi- 
madas para u pode ser da maior utilidade por permitir simplificar em considerável extensão aquele 
problema analítico cujas soluções seriam de outro modo praticamente inviáveis. 


APÊNDICE 


Definindo o coeficiente de resistência do escoamento no modo habitual 


f=—-B>—— (13) 


AA mm 


calcula-se imediatamente o coeficiente de resistência para um canal de secção rectangular e qualquer 
a : a 
F = » a partir das expressões apresentadas para Um. 


Em função do número de Reynolds à expressão de f é como facilmente se verifica dada por 


ae 32 
k R (1 + ) 

q cd 

a b b «A . , . Um 
Sendo D=4—— =4 -—— oq diâmetro hidráulico e Kk = — 

a+b dd + a =. 

F u 
f E 1 b ; 
Apresentam-se seguidamente os valores de —— em função de a =— (b<a) 
Te a 


(13) Em escoamento uniforme, num canal das condições do modelo, se P; e P> forem as pressões na secção 1 e 2 
respectivamente, à distância L e as cotas h; e hz é ôbviamente 
P, — 


ns P, 


(hy—ho) 1 
L L 


TECNICA N.º 362 71 


Coeficiente de resistência f para escoamento 


uniforme em canais de 


laminar 


4 
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tratada com permutadores de iões Lewatit. 


e 
porra 


Água como esta, muito rica em ácido humínico, também pode ser utilizada para alimentar caldeiras, se for 


E preciso ir buscar aqui a água para as suas instalações de caldeiras? 


Pode fazê-lo sem receio. Com os permu- 
tadores de iões Lewatit é possivel trans- 
formar água superficial, suja, contendo 
ácidos humínicos e detergentes, numa 
excelente água para alimentação de cal- 
deiras 

A Bayer desenvolveu numerosos tipos de 
PLewatit para descalcificar, desmineralizar 
e dessilicificar a sua água. Especialmente 
as resinas macroporosas de estrutura 
esponjosa apresentam grande resistência 
à contaminação irreversivel por substân- 
cias orgânicas. Os nossos tipos MP,macro- 
porosos, são capazes de absorver ácidos 
humínicos e detergentes no decorrer dum 
ciclo e de os ceder novamente só na rege- 
neração. A estrutura esponjosa dos nossos 
tipos macroporosos MP pode ser reco- 
nhecida, na prática, pelo seu excelente 
comportamente na presença de águas 
contaminadas por substâncias orgânicas. 


Sob o microscópio electrónico torna-se 


visivel a razão deste vantajoso comporta- 


mento: 


Dos poros principais com um diâmetro 
de alguns milhares de Angstróm deriva 


um sistema finamente ramificado de ca- 


nais. Os macroporos conduzem assim 
a uma superficie interior de fácil acesso. 
isto proporciona uma rápida difusão das 
partículas carregadas (iões) e das não 
carregadas (ácidos huminicos e detergen- 
tes), dissolvidos na água, nos centros de 
actividade permutante da resina, 


Este mesmo sistema heterogénio de poros 
é condição essencial para uma eluição 
irrepreensível das substâncias orgânicas 
combinadas, no processo de regeneração. 
A depuração duma resina macroporosa, 
das substâncias orgânicas absorvidas, 
efectua-se no decorrer duma regeneração 


normal, com soluto alcalino. Os permuta- 
dores de iões macroporosos não necessi- 
tam por isso normalmente de tratamento 
especial. 

A Bayer fornece, além dos permutadores 
de catiões e de aniões, ainda resinas 
absorventes. 

Tem um problema técnico? Confie na 
experiência da Bayer! 


Bayer Alemanha 


PSY lit, 


Representante em Portugal: 

S. A. R. L..Quimicor 

Lisboa, Rua Sociedade Farmacêutica, 3, 
Tel. 42194, 

Porto, Rua de Santos Pousdda, 441, 

Tel. 54141 
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O mesmo 
que se obtém com 
um telefone normal 
através de um 
circuito físico, pode ser 
agora obtido com um 
feixe hertziano tipo 

ATE (700) transistorizado 
sem usar qualquer linha física 


Quer seja de telefone para telefone 
ou telefone para central ou ainda central 
para central, o equipamento ATE (700) 

é o mais recente e de melhor utilização 
nos centros rurais. Melhor pelo seu baixo 
consumo. pois consome apenas uma 

carga de pilhas por ano. Melhor porque 
sendo um rádio telefone tem as mesmas 
possibilidades que tem um telefone, 
inclusivé facilidades de funcionamento 
com linhas partilhadas num circuito 

de transmissão. 

Melhor ainda porque não tem válvulas — é 
completamente transistorizado, de modo 

a ter a maior duração possivel. 

É de instalação fácil, completamente 
estanque podendo trabalhar entre 
temperaturas de — 20" C a + 55ºC. 
Fabrica-se nas bandas alta e baixa de 

VHF (54/88 Mc/s). (150/184 Mc/s) 

te sdo ainda a possibilidade de ser usado 
com um amplificador de Rádio Frequência 
para um aumento do raio de propagação. 
Estes são portanto os motivos que nos 
levam a aconselhar a adopção de tal 
equipamento em áreas de dificil acesso 
para manutenção. 

Para mais detalhes os Serviços Técnicos 

da AEP estão à inteira disposição 

de V. Ex.**. 
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em, 


Tensão de segurança 3. 000 Kgf / cm2 
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AÇOS DE ALTA RESISTÊNCIA PARA BETÃO ARMADO 


que pela eficiência técnica e económica a que conduzem, são 
um indispensável clemento na moderna construção civil 


Para qualquer obra em betão armado cossuite: 
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A INVESTIGAÇÃO DA ESTRUTURA DAS SOLUÇÕES DE ELECTRÓLITOS 
POR MEIO DA DIFRACÇÃO DOS RAIOS X 


RESUMO 


Indica-se o objectivo actual da instalação para a difrae- 
ção dos raios X recentemente montada no Laboratório de 
Física do Instituto Superior Técnico * a qual se destina 
ao estudo de estruturas moleculares, Referem-se algumas 
circunstâncias do seu projecto e faz-se a descrição dos ele- 
mentos da montagem nos sens aspectos essenciais. 


INTRODUÇÃO 


São consideradas aqui apenas as soluções aquo- 
sas dos sais inorgânicos de catiões metálicos 
com a estrutura electrónica dos gases raros; 
serão postas de parte as interacções entre o 
ião H* e as moléculas do solvente e não será 
analisada a interacção entre o anião e as molé- 
culas de água. 

A ideia de que os iões em solução coordenam 
em torno de si as moléculas de água do solvente 
tem sido utilizada, pelo menos qualitativamente, 
na descrição de um grande número de proprie- 
dades das soluções electrolíticas, mas o desco- 
nhecimento da forma como as moléculas de água 
estão associadas ao ião central leva por vezes à 
utilização de esquemas demasiado simplistas: os 
cálculos baseados nesses mecanismos conduzem 
para uma mesma solução electrolítica a um nú- 
mero de moléculas do solvente associadas aos 
iões que é variável conforme a propriedade física 
ou química que se pretende explicar. 

A reorganização da estrutura da água devida 
à presença de iões dissolvidos pode ser posta em 
evidência pela alteração do volume da solução 
aquosa relativa a cada molécula grama de elec- 
trólito que se dissolve a mais (1), pela alteração 


por MANUEL ALVES MARQUES 
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Bolseira do Instituto de Alta Cultura 


SYNOPSIS 


An X-ray diffraction arrangement for the study of 
molecular structures is describea in its main points, 
Some circunstances of the project are pointed out. 


nas mesmas condições da compressibilidade (1), 
(2), da capacidade calorífica (3), da constante 
dieléctrica (4), do índice de refracção (5), da vis- 
cosidade (1), da condutividade eléctrica (6), da 
solução aquosa, pela alteração da frequência de 
ressonância magnética dos núcleos do hidrogé- 
nio ou dos núcleos de oxigénio (7) (8) das molé- 
culas de água, pela medição do calor e da entro- 
pia (10) de dissolução, etc. 

O confronto do número de hidratação que é 
o número de moléculas de água ligadas ao ca- 
tião determinado por alguns métodos experi- 
mentais — não se pretendem mencionar aqui todos 
os numerosos trabalhos referentes à hidratação 
dos iões em solução até agora publicados — ilus- 
tra a discordância referida atrás. 


Números de Hidratação 


MÉTODOS EXPERIMENTAIS Lit Mg*t K+ 
Medição da capacidade ca- 

lotitica 13) sw é sed 19 E DR À 
da velocidade de propa- 

gação das ondas hiper- 

sônicas (11). ... cv c. 10a 6 4 15a12 


da constante dieléctrica (4) 6 5 cm 


* Alguma aparelhagem foi adquirida com uma verba concedida pelo Instituto de Alta Cultura mediante um 


acordo estabelecido entre esta entidade e a Shell Portuguesa, 
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E que a introdução dos iões na água deter- 
mina nalguns casos dois efeitos contrários e às 
vezes não é fácil destrinçar a influência de cada 
um deles: um, ao nível da primeira camada de 
hidratação, promove uma solidificação do líjuido 
fixando ao catião as moléculas de água dessa 
camada, o outro, ao nível das segunda, terceira 
camadas rompe a estrutura tetraédica da água 
líquida e promove assim um aumento de flui- 
dez (12). Estes dois efeitos influem nos re- 
sultados que os diversos métodos experimen- 
tais permitem obter de maneira muito dife- 
rente. 

Como as moléculas de água directamente liga- 
das aos catiões metálicos constituem um edifício 
suficientemente estável para se poderem definir 
estados de vibração própria a espectrografia do 
efeito de Raman — e eventualmente a do infra- 
vermelho — podem fornecer importantes resul- 
tados relativos às forças que constituem o com- 
plexo, às suas propriedades de simetria e à sua 
vida média. Estes complexos constituídos por 
uma camada de moléculas de água aderentes ao 
catião metálico foram estudados pelo efeito de 
Raman pela primeira vez por A. da Silveira (13) 
em 1932 o qual conseguiu observar as ban- 
das que individualizam o complexo octaédrico 
Me*(OH:), de moléculas de água nas soluções 
de sais de alumínio de magnésio de zinco e de 
cobre. Mais tarde, em 1961 (14), estabele- 
ceu definitivamente a estrutura dos complexos 
Mg”*(OHs) e ALt(OH>) e em 1965 (15) pu- 
blicou investigações sobre as soluções de sais 
de berílio que revelam as bandas de Raman 
devidas às vibrações dos complexos Be**(OHy),. 
Investigações mais recentes (16) permitem admi- 
tir a coexistência de complexos Be?*(OH:) e 
de complexos Be?+ (OHs);*+ em solução. Estes 
resultados experimentais foram confirmados pri- 
meiro por J. P. Mathieu (17) que estudou pelo 
efeito de Raman os cristais e as soluções aquo- 
sas de certos sais; e mais recentemente por R. 
Lafont (18), R. Hester e R. Plane (19), e por D. 
E. Irish (20). 

É a observação de um espectro de Raman 
idêntico nas soluções de diversos sais de um 
mesmo catião, espectro de Raman que é dife- 
rente de catião para catião, que permite afirmar 
a existência de uma estrutura intermolecular es- 
tabelecida por esse catião. A persistência das 
bandas de Raman quando a concentração dimi- 
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nui garante a presença da estrutura intermole- 
cular nas soluções diluídas e simultâneamente 
afasta o hipótese de se tratar de um fenómeno 
característico das altas concentrações porque a 
intensidade das bandas diminui sensivelmente 
segundo a mesma lei que a concentração do 
catião. Como os cristais hidratados se depositam 
da solução saturada decerto que a configuração 
interna da solução na vizinhança de uma grande 
parte dos catiões presentes consiste numa ca- 
mada de moléculas de água ligadas a estes 
iões tal como nos cristais hidratados. A di- 
fracção dos raios X pelos sais hidratados (21) 
dos catiões respectivos fornece pois um ponto 
de partida para o número de moléculas de água 
adjacentes ao catião, a cuja aderência a obser- 
vação pela espectrografia do efeito de Raman 
permitiu atribuir uma vida média relativamente 
longa, e fornece a distância entre o catião e as 
moléculas de água do mesmo complexo. Estas 
distâncias — raios de equilíbrio dos complexos — 
são ligeiramente discordantes nos diversos sais 
cristalizados do mesmo catião porque no estado 
sólido a cooperação entre os diversos consti- 
tuintes da rede cristalina é mais efectiva e per. 
mite assim que os aniões deformem ligeiramente 
os complexos. 

São pouco numerosas as determinações feitas 
até hoje sobre as soluções electrolíticas com o 
objectivo de determinar as distâncias entre o 
catião e as moléculas de água (22) (23). Ne- 
nhumas incidiram sobre uma estrutura firme- 
mente estabelecida previamente pela espectro- 
grafia do efeito de Raman. A montagem experi- 
mental que se descreve a seguir permitirá obser- 
var em solução por meio da difracção dos raios 
X os complexos da forma Me(OH:) cuja indi- 
vidualidade foi investigada já por métodos es- 
pectográficos. 


MONTAGEM EXPERIMENTAL 


É conhecido que as figuras obtidas por meio 
da difracção dos raios X pela matéria no estado 
liquido são diferentes das figuras de difrac- 
ção produzidas pelos gases e diferentes das 
figuras de difracção produzidas por um cristal: 
uma chapa fotográfica em que incidiu um feixe 
de radiamento X monocromático difractado por 
uma fina camada de um gás monoatómico mos- 
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tra uma constante diminuição da densidade de 
enegrecimento a partir do centro onde incidiu o 
feixe directo até à periferia. Na figura da difrac- 
ção de um feixe de raios X por um cristal pul- 
verizado observa-se, a partir de uma certa dis- 
tância do centro onde incidiu o feixe directo 
uma sucessão de anéis concêntricos cujas inten- 
sidades variam por variar a concordância de fase 
dos raios X que interferem em cada ponto do 
alvo. A estrutura electrónica do cristal permite 
calcular satisfatóriamente as intensidades obtidas 
sobre a chapa fotográfica por meio da teoria 
ondulatória dos raios X. 

Se o meio material que produz a difracção 
dos raios X for um líquido observa-se natural- 
mente uma figura de simetria circular que não 
coincide, apesar de se notarem algumas seme- 
lhanças, com a série de anéis de difracção pro- 
duzidos pelo respectivo cristal pulverizado. Os 
anéis, em menor número, são mais difusos e 
vão-se esbatendo até desaparecerem à medida 
que o seu raio aumenta. Nas diferentes formas 
de considerar as posições relativas médias dos 
átomos no líquido, formas que permitam explicar 
satisfatóriamente os diagramas de difracção, con- 
sistem os diversos modelos do estado líquido da 
matéria. Estes modelos conservam vestígios mais 
ou menos pronunciados de uma estrutura crista- 
lina cuja extensão e subsistência têm sido larga- 
mente discutidas nos últimos anos. 

É possível arquitectar modelos intermolecula- 
res de um determinado líquido bastante dife- 
rentes entre si que, no entanto conduzam por 
meio da análise de Fourier a diagramas de difrac- 
ção praticamente indistinguíveis dentro dos limi- 
tes da precisão experimental. Apesar desta lati- 
tude da indeterminação do modelo do líquido a 
partir dos resultadcs experimentais fornecidos 
pela difracção dos raios X, tem sido estudada 
por este método uma extensa série de substân- 
cias no estado líquido (24) (25) (26) — os gases 
raros, os metais, os compostos orgânicos, a água, 
os sais, as soluções aquosas e não aquosas, as 
substâncias vitreas. A estrutura de muitos destes 
liquidos é quase estabelecida à priori; reser- 
vam-se as experiências de difracção para deter- 
minar a grandeza das distâncias intermoleculares, 
ou para decidir, entre um pequeno número de 
arranjos intermoleculares cujas figuras de difrac- 
ção se prevêem bem diferentes, qual é a forma 
interna do líquido investigado. 


TÉCNICA N. 362 


Descrição 


O feixe de raios X é estritamente monocro- 
mático para permitir a interpretação pormeno- 
rizada da figura de difracção obtida no alvo. 
Usa-se como monocromatizador a lâmina encur- 
vada de um cristal. 

A montagem resulta geralmente mais mane- 
jável quando a difracção se executa em planos 
horizontais. A superfície livre do especímen será 
por isso vertical. No estudo da água e das so- 
luções aquosas têm sido usadas até hoje quase 
exclusivamente dois tipos de montagem do espe- 
cimen: a) montagem de jacto livre vertical em que 
o líquido é exposto directamente ao feixe de raios 
X sem a interposição da parede da célula que o 
contém (23) (27); b) montagem em que o Ji- 
quido está contido no interior de uma célula de 
paredes transparentes aos raios X (22) (28) (29). 
O primeiro tipo tem a vantagem evidente da 
pureza da amostra — não é necessário descontar 
a difracção devida à parede da célula. O segundo 
tipo tem a vantagem de permitir uma regulação 
cómoda da expessura da solução a usar e de per- 
mitir expor ao feixe de raios X uma região ex- 
tensa: um jacto livre dificilmente se regula até um 
diâmetro de 0,4 mm — quando se tenta obter um 
diâmetro menor que este começa a dividir-se em 
gotinhas —e para aproveitar todo o feixe de 
raios X teria de se utilizar uma cortina de peque- 
ninos jactos o que não é cómodo. Optou-se pelo 
segundo tipo — célula plano paralela (fig. 4) — 
exigindo-se então um ensaio a branco, com a 
célula vazia para corrigir a figura de difracção 
obtida. Esta montagem foi já usada por F. De- 
bot (28) para estudar os complexos Cl;Zn?— nas 
soluções aquosas de Cl:Zn. As paredes da célula 
eram películas de nitrocelulose com a espessura 
de menos de 0,1 mm. G. Brady (29) usou uma 
célula com parede de berílio com 0,125 mm de 
espessura para estudar as soluções electrolíticas. 

A utilização de contadores Geiger-Miiller na 
exploração do alvo onde se forma a figura de 
difracção tem a grande vantagem de rejeitar o 
radiamento de fluorescência do especimen, cuja 
frequência é mais baixa que a frequência do 
radiamento X do feixe incidente (29): os im- 
pulsos correspondentes à região espectral da 
fluorescência do especimen são muito mais 
fracos e não atingem por isso o limiar do regis- 
tador de impulsos que só terá em conta os im- 
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